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An analytical expression is developed 
for describing the thermotropic behav- 
ior of membrane bilayers as studied by 
Raman spectroscopy. The expression 
is derived from a two-state model of 
the main gel to liquid crystalline phase 
transition in lipid bilayers. Experimen- 
tal data for a variety of diacylphos- 
phatidylcholines and their derivatives 
have been fit by least squares to the 
two-state expression to within cur- 
rently achievable measurement error. 
Numerical techniques have been de- 
veloped for placing bounds on the 



parameters of the two-state model in 
situations of sparse data in the phase 
transition region. By fitting the model 
to the measured spectroscopic data, es- 
timates of the extent of cooperativity 
in the phase transition can be obtained 
in a systematic manner. 

Key words: logistic function; mem- 
brane; phase transition; Raman spec- 
troscopy; two-state model. 

Accepted: December 3, 1986 



Introduction 

In recent years the use of Raman spectroscopy 
for the study of the phase transition and the order/ 
disorder behavior of both intact and reconstituted 
membrane bilayers has gained dramatically in pop- 
ularity through the wide range of applications 
amenable to spectroscopic techniques [1]\ Al- 
though calorimetric measurements can provide ac- 
curate information on enthalpy changes that occur 
during phase transitions, they provide little insight 
into the nature of the molecular reorganization ac- 
companying the phase transition. Since Raman 
spectra are particularly sensitive to the conforma- 
tional, packing, and dynamical changes involving 
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the hydrocarbon chains, bilayer reorganizations 
within, for example, the hydrophobic region of the 
membrane, may be monitored directly as a function 
of temperature. These observations can be made in 
an aqueous medium, without either the introduc- 
tion of bulky molecular probes or elaborate sample 
preparation. By using low laser excitation powers 
(15-200 mW) at the sample, external perturbations 
to the system are minimized. Because of the sensi- 
tivity of the spectra to conformational changes, Ra- 
man spectroscopy has also been useful in 
elucidating the effects of intrinsic and peripheral 
components such as sterols and proteins on the 
thermal properties of membranes. In a few systems, 
reconstituted, naturally-occurring membranes have 
been studied and have provided useful information 
obtainable only from spectroscopic examination 

[2]. 

Although the vibrational assignments for the 
spectra of lipid bilayers are complex, they are un- 
derstood sufficiently well to allow discussion of 
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membrane reorganizations in terms of the observed 
spectroscopic frequencies, intensities and line- 
widths of a variety of spectral features reflecting 
the headgroup, interface, and acyl chain regions of 
the lipid components. A particularly useful ap- 
proach for presenting the thermal data has been in 
terms of spectral peak height intensity ratios deter- 
mined from either the 2900 cm -1 carbon-hydrogen 
or 1100 cm -1 carbon -carbon stretching mode re- 
gions. In general, these data are primarily sensitive 
either to alterations in the lateral interactions or to 
intermolecular order of the acyl chains and intra- 
chain trans -gauche conformational changes, respec- 
tively. 

In this paper we propose an analytical expres- 
sion, based upon a two-state model of the primary 
phase transition phenomenon, which is useful for 
systematically describing the form of the tempera- 
ture profiles. We also describe an approach, based 
on a nonlinear, least squares analysis, for extracting 
relevant thermodynamic information from the tem- 
perature dependence of the observed spectra. This 
approach is robust in handling the range of data 
encountered in studies of membrane phase transi- 
tions. 



Thermodynamic Model 

Within the last 10 years a number of models have 
appeared describing within either a thermody- 
namic or statistical mechanical framework various 
features of the main phase transition in membrane 
bilayers. Most of these models have been critically 
reviewed by Nagle [3]. An additional review by 
Israelachvili, et al. [4] discusses the physical princi- 
ples of membrane organization and includes de- 
scriptions of membrane-protein interactions and 
nonbilayer lipid structures. We concern ourselves 
here with the simplest of these phase transition 
models; namely, that of thermodynamic equi- 
librium between noninteracting domains in either 
one of two states [5]. We have been guided by the 
philosophy that the accuracy with which a simple 
model can quantitatively describe the experimental 
data serves as a guide to assessing the required fea- 
tures of a more complex model, to comparing and 
evaluating two or more models, and to defining the 
accuracy requirements of additional measurements. 

The derivation of the thermodynamic model is 
presented in Appendix A in detail in order to dis- 
tinguish clearly between numbers of molecules and 
numbers of domains in each of the two states since 
the observed spectral intensities are determined by 
the number of molecules in each of the states while 



the temperature dependence of the spectral intensi- 
ties is determined by the distribution of domain 
sizes. The mathematical development follows 
closely the description by Hill of the two-state 
model for small systems [6]. The resulting analyti- 
cal expression for the temperature dependence of 
the Raman spectra are identical to those published 
by Dluhy, et al. [7] using the Zimm-Bragg theory 
originally developed for describing polypeptide de- 
naturation [8]. The difference between the thermo- 
dynamic and the Zimm-Bragg developments is in 
the interpretation of the parameters derived from 
the analysis. We prefer the thermodynamic devel- 
opment because of the simplicity of its physical in- 
terpretation. For so long as the data are described 
by a simple model, it seems little more than a mat- 
ter of taste as to how the parameters are to be inter- 
preted. Only when the data can be seen to contain 
additional information would it be fruitful to apply 
a more complex model. Finally, though the deriva- 
tion is applied in this paper to Raman spectroscopic 
data, it is generally applicable to all spectroscopic 
studies of such cooperative phenomena. 

In the two-state model, a single bilayer is consid- 
ered to consist of domains of varying numbers of 
lipid molecules n with each domain existing in one 
of two states labeled A and B. The distribution of 
domain sizes and the requirement that all molecules 
within a domain exist either in state A (the low- 
temperature gel state) or in state B (the high-tem- 
perature liquid crystalline state) determine the 
cooperativity exhibited by the bilayers during the 
phase transition between the two states. The distri- 
bution of domain sizes is assumed to be tempera- 
ture independent, while the number of domains in 
either state is temperature dependent. Although a 
phase transition in the strictest sense is excluded for 
a distribution of domains of finite size, the discon- 
tinuous nature of the phase transition will be ap- 
proximated for large values of n . 

The temperature dependence of the intensity Iv 
of a spectral feature (Raman or otherwise) that de- 
pends upon the distribution between two states of 
the molecules being observed is given by eq (A32) 
in Appendix A: 



■Mi; 



:+ 



(1+eTO+e-') 



(1) 



The dimensionless, reduced temperature t is de- 
fined as t = (T—T )/D, where T is the midpoint of 
the phase transition and the scale factor D is related 
to the enthalpy of transition per mole of lipid 
molecules, A//, and the effective domain size, n efr , 
by: 
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D = {Rlt)/bHna 



(2) 



As discussed in Appendix A, eq (1) represents an 
approximation to a summation over a distribution 
of domain sizes. To test the accuracy of this ap- 
proximation, sets of hypothetical data were created 
using a lognormal distribution [9] of domain sizes 
from which Raman intensities were calculated for 
values of N A and N B , the numbers of molecules in 
states A and B respectively. The values of N A and 
N B were calculated for the lognormal distribution 
using eqs (A21) and (A24). The lognormal distribu- 
tion was chosen because it excludes domains of 
negative size and allows for a pronounced skewing 
to high values of n . This skewing causes various 
measures of the center of the distribution to differ 
significantly and provides a sensitive test of a sin- 
gle-term approximation to a sum over a distribu- 
tion. The data sets corresponded to distributions of 
varying breadth which, when fit by the method of 
least squares to eq (1), gave the same value of n efr . 

In table 1, five measures of the center of a distri- 
bution: the mode, M, which is the most populous 
value of n, the mean value of n, <n >, the root 
mean square of n, V</i 2 >, the median, n(50%), 
and the value of <# 2 >/<«> are compared with 
the value of /i e fr returned by the least squares fit for 
three distributions of varying spread. It can be seen 
that as the spread of the distribution narrows, the 
agreement of all of these measures improves. Even 

Table 1. Comparison of a single term fit to three distributions. 



though, for the broadest distribution, the values of 
these measures covered more than an order of 
magnitude, eq (1) gave a satisfactory fit with a rela- 
tive standard deviation of 2%. 

Other conclusions may be drawn from the data 
presented in table 1. The value of « efT obtained from 
the least squares fit was best approximated by the 
root mean square value V<« 2 >, and the median, 
«(50%), rather than by <« 2 >/<«> as we had 
initially expected. This is a function, however, of 
the density of data in the vicinity of the transition 
interval. If the data range is narrowed to include 
only the transition interval and if the number of 
data in that region is increased, then the limit of 
<n 2 >/<n> is approached in agreement with the 
observation in Appendix A that this measure gives 
the correct temperature dependence at the asymp- 
totes and at the inflection point when used in a 
single term approximation to the sum over domain 
sizes. Another pronounced effect of the spread of 
the distribution was to give a more gradual conver- 
gence to the asymptotic limits and thus, the appear- 
ance of sloping rather than horizontal asymptotes. 
This can be seen in figure 1 in which the hypothet- 
ical intensity ratio profiles are drawn for each of 
the three distributions listed in table 1. The sloping 
asymptotes arise from contributions from domains 
of small size as originally suggested by Sturtevant 
[10]. 

A somewhat better fit to broad distributions can 
be obtained by modifying eq (1) in order to allow 
for nonhorizontal asymptotes, namely: 



Mode 



1 30 90 



Range in n* 1-11,679 1-3,963 21-387 

Distribution <n> 36.3 79.2 98.2 

Parameters 

V<« 2 > 113.7 109.5 101.0 

<n 2 >/<n> 356.3 151.3 104.0 

median 117 110 101 



Least Squares Fit n d 
to Eq 1 



100 100 100 

2.1% 1.2% 0.2% 



toEq3 «eff 136.2 113 101 

s 1.0% 0.8% 0.2% 

slope c 0.46% 0.24% 0.02% 

"Range: values of n for which the population >10 -8 x population of 
mode. 

b s: relative (to the value of A-B) standard deviation of the least squares 
fit. 

c slope: value of baseline slope relative to the value of A-B. 



_ I Av +I av (T-T ) L 
IV ~ (1+eO + 



(1 + e-') 



(3) 



From table 1 it can be seen that inclusion of a 
nonhorizontal baseline in the fit of the hypothetical 
data improved the quality of the fit, as judged by 
the standard deviation, and increased the estimate 
of the value of n slightly, though not nearly enough 
to bring it to its expected value of <« 2 >/<n>. 

In the absence of an analytical expression such as 
eq (3), it has been customary to estimate T and the 
van't Hoff enthalpies from observed temperature 
profiles using straight line approximations to the 
three distinct regions of the profiles. Thus, AT has 
been defined as the temperature interval between 
the intersection of a line drawn tangent to the 
phase transition curve at its inflection point and the 
two lines drawn tangent to the high and low tem- 
perature asymptotes. This definition of A T is illus- 
trated in figure 2. The relationship of A T to the 
scale parameter D of eq (3) can be calculated di- 
rectly to give AT=4Z), ignoring terms of order 
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Figure 1-A graphical comparison of 
the temperature dependence of hy- 
pothetical Raman intensity ratios 
for three domain size distributions 
with modes M=\, 30, and 90 and 
n c{r = 100. 



TEMPERATURE 
(ARBITRARY UNITS) 



Neft = 4kTo 2 
AH AT 




Figure 2-Relationship between the 
parameters of eq (5) and the graphi- 
cal analysis of phase transition 
curves. 



TEMPERATURE 
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(Ibv—IdvVilBv—lAv)- From the definition of D given 
by eq (2), 



Ar=(4^r 2 )/A^n efT . 



(4) 



Comparison of eq (4) with similar expressions 
given by Sturtevant [11] indicates that AT" corre- 
sponds to the AT used for calorimetric determina- 
tions of van't Hoff enthalpies. Care should be taken 
in the interpretation of n ef{ , however, since esti- 
mates of the slope made as close to the inflection 
point as possible will give values of « eff approach- 
ing <« 2 >/<«> which, for broad distributions of 
domain sizes, can be an order of magnitude higher 
than </i >. 

In the Raman spectroscopic experiment tempera- 
ture profiles are collected over periods of several 
hours. Since fluctuations in the intensity of the inci- 
dent laser radiation lead to corresponding fluctua- 
tions in the observed Raman intensities, the ratio of 
the intensities of two spectral features at different 
scattering frequencies provides a more precisely 
defined temperature profile which is relatively in- 
dependent of the experimental environment. More- 
over, if the vibrational features are chosen to be 
sensitive to the ordering of the lipid hydrocarbon 
chains, the temperature dependence of the intensity 
ratio will be more pronounced than for the individ- 
ual intensities themselves. 

Forming the ratio r i2 =I\/I 2 and substituting eq 
(3) for I x and I 2 , we obtain after rearrangement, 



, A+a(T-T )B+b(T-T ) 
rn ~ (1+eO + (l + e-0 

where, to first order in 1^,/Iav or I^/Ibv, 

-M2 -M2\ \ l Al) J J B2 



(5) 



and Tl = T-Dkife\ 



(T-TS) 
D ' 



(6) 



In the expression for To, the contributions from 
hi and I al have been ignored. A comparison of eqs 
(3) and (5) shows that the intensity ratio has ap- 
proximately the same functional dependence on 
temperature as the individual intensities them- 
selves, except that To, the apparent midpoint of the 
phase transition, has been shifted from the true 
value by an amount Dlnilsi/lAi)- If» instead of r n , 
the temperature behavior ofr 2l =I 2 /I\ is analyzed, a 



new apparent midpoint, To=T — DXn(I Bl /I Ai ), is 
obtained and a comparison of the temperature de- 
pendence of r 2 \ with r n shows an apparent shift of 
the transition temperature of D\n(B/A). 

Eqs (3) and (5) are the expressions we have used 
for the analysis of the Raman intensity profiles. 
Spectral data are fit by the method of least squares 
to Raman intensities or intensity ratios. From these 
fits, estimates of transition temperatures and van't 
Hoff enthalpies are obtained in a systematic man- 
ner. For some systems in which this analysis has 
been applied, the sharpness of the phase transition 
requires a particularly robust least squares calcula- 
tion. Therefore, the details of the least squares anal- 
ysis are presented in Appendix B. 



Analysis of Temperature Profiles 
Derived From Raman Intensity Ratios 

The analytical approach outlined in the preced- 
ing section and in the appendices to this paper have 
been applied to dipalmitoylphosphatidylcholine 
(DPPC) as well as to several cyclopentanoid ana- 
logues of DPPC [12]. In the following discussion, 
the results of some of these analyses will be de- 
scribed for the purpose of examining the limitations 
of the two-state model and the ways in which these 
limitations are made manifest by the statistical anal- 
ysis. 

Although several Raman spectral regions are 
suitable for deriving temperature profiles from ap- 
propriate intensity ratios to assess bilayer reorgani- 
zations, we will use the 2900 cm -1 C-H stretching 
mode region for applying the present model. This 
congested and complex spectral interval, shown in 
figure 3, for the gel and liquid crystalline states of 
DPPC has been discussed numerous times (see, for 
example, ref. [1] and references contained therein). 
Temperature profiles can be conveniently con- 
structed from changes in the peak height intensity 
ratios for the 2850, 2885 and 2940 cm" 1 feature. 
These spectral transitions are assigned, respec- 
tively, to the acyl chain methylene C-H symmetric 
stretching modes, the methylene C-H asymmetric 
stretching modes and, in part, a Fermi resonance 
component of the chain terminal methyl symmetric 
C-H stretching mode [1]. 

Figure 4 illustrates the analysis of a typical phase 
transition profile for multilamellar dispersions of 
DPPC using the spectral intensity ratios con- 
structed from the C-H stretching mode region. The 
values of A, B, a and b listed in the figure are 
the values of the parameters of eq (5) along with 
their estimated uncertainties. Only one datum, at 
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Figure 3-Raman spectra of dipalmitoylphosphatidylcholine in 
(A) the gel state at 35.6 *C and (B) the liquid crystalline state 
at 42.3 *C. 



40.67 °C, fell within the transition interval, and 
consequently A T could not be determined. Illustra- 
tive of the robustness of the analysis, it was as- 
signed a maximum value based on the temperature 
separation between the data on either side of 
40.67 °C (39.58 and 41.78 °C) and on the standard 
deviation of the fit using eq (Bll) in Appendix B. 
The transition data range given in figure 4, 1.10 K, 
is the value of the difference between the measured 
temperatures below and above T and should be 
viewed as spanning the uncertainty in T . Ninety- 
five percent confidence limits for T can be ob- 
tained from the standard deviation of 7^ returned 
by the least squares fit, 0.08°, by multiplying by 
12.7, the value of the 95% confidence limits of a 
student's t distribution with one degree of freedom, 
giving 1.02 K in agreement with the transition data 
range. As might be expected, the value of T as 
determined by the least squares fit is quite sensitive 
to the value estimated for A T. The estimated in- 



verse shift presented in figure 4 is calculated from 
Dln(B/A) and is the shift in T which would have 
been obtained if 7(2885)/7(2940) had been analyzed 
instead of 7(2940)//(2885). The vertical bars at 
each data point in figure 4 are twice the value of 
the uncertainties in the intensity ratios as estimated 
from the standard deviation of the fit by eq (Bl). 
The three curves passing through the data values 
correspond to the profile calculated from eq (5) 
along with plus or minus two standard deviations. 
Along the horizontal axis of figure 4 is the two 
standard deviation curve alone. The peak in the 
uncertainty curve at T Q signifies that T Q is being 
determined almost exclusively by the value of the 
intensity ratio at the one datum in the transition 
interval relative to the separation between A and 
B. 

Figure 5 illustrates several of the limitations of 
the simple two-state model. Four phase transition 
curves for DPPC from our laboratory using data 
from different investigators, different spectrome- 
ters, and over a two-year interval were compared. 
Even though the data were separated by 2° inter- 
vals, and only one, but always one, datum fell 
within the transition interval, the analysis provided 
a surprisingly reproducible phase transition tem- 
perature of 40.38 to 41.09 °C, as compared with the 
more accurate value of 41.55 °C measured calori- 
metrically by Albon and Sturtevant [5]. The maxi- 
mum value of AT varied from 1.42 to 1.91 °C and 
reflected varying temperature data intervals and in- 
tensity measurement precision. 

Using the calorimetric data for DPPC measured 
by Albon and Sturtevant, namely, A//=8.50 kcal/ 
mole and r =314.70 K, and the value obtained for 
A7\ minimum values of the cooperative unit of 48 
to 65 were calculated from eq (4). If values for the 
cooperative unit of 100 to 1000, which have been 
estimated for the phase transition in DPPC, are to 
be obtained from the analysis of Raman intensity 
profiles, measurements must be made with temper- 
ature increments of 0.1 °C in the phase transition 
region and with temperature control on the order 
of 0.01 °C. 

If, for the moment, the data for the second run of 
sample # IB are excluded from figure 5, the values 
of A, B, a and b all agree to within three standard 
deviations as estimated by the least squares fit. For 
the temperature profile of sample #1B, which was 
measured within 24 hours of the profile of sample 
#1A, the values of A , B, a and b all varied by more 
than five standard deviations from the average of 
the other three profiles. We interpret this observa- 
tion as an indication that the sample had not yet 
annealed into its preferred low temperature state. 
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Figure 4-Temperature profile for the ratio of the intensity of the Raman 
spectral feature associated with, in part, the Fermi resonance component 
of the acyl chain terminal methyl symmetric C-H stretching mode at 
=2940 cm -1 and the methylene asymmetric C-H stretching modes at 
=2885 cm -1 of an aqueous dispersion dipalmitoylphosphatidylcholine 
(DPPC). The central curve represents the least squares fit of eq (5) to the 
spectral intensity data. The upper curve represents the calculated curve 

plus two standard deviations and 
the lower curve minus two stan- 
dard deviations. The curve drawn 
along the bottom axis is the two 
standard deviation curve alone. 
The vertical bars are estimates of 
the uncertainties in the experi- 
mental ratios and are based on the 
weights of the ratios and the 
overall standard deviation of the 
fit. The values of the parameters 
obtained from the fit are as fol- 
lows: A =0.474±0.006, 5=0.740 
±0.012, a =0.0074±0.OO08 and 
b = 0.0036 ± 0.0023, T = 40.38 
±0.08 "C and A7'<1.76 , C. The 
transition data range = 1.10 *C. 
The estimated inverse shift 
= — 0.10'C. The standard devia- 
tion 5=1.70% relative to the 
value of r at the midpoint of the 
transition interval. 
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Figure 5-A comparison of four tem- 
perature profiles for dipalmi- 
toylphosphatidylcholine (DPPC). 
The spectra were measured over a 
period of two years by three differ- 
ent investigators. Spectra for Sam- 
ples No. 1 and No. 3 were recorded 
with a scanning spectrophotometer 
whereas spectra for Sample No. 2 
were determined with a spec- 
trograph. All three samples were 
prepared by the individual investi- 
gators using different sources for 
the lipid. The range in T was 
40.4 'C to 41.1 'C. 
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Its bilayer structure appeared, on comparison with 
the values of the gel state intensity ratios for the 
other samples, to be as ordered at low tempera- 
tures, although it did not exhibit the gradual in- 
creasing disorder associated with the positive value 
of a. Its phase transition temperature, though 
somewhat higher than the others, did not change 
significantly nor did the transition interval to 
within measurement precision. The transition con- 
cluded at a much less disordered state than the 
other samples and then, corresponding to a large 
value of b, continued to disorder at a much greater 
rate than the others until reaching a comparable 
intensity ratio approximately 10 °C above the phase 
transition temperature. Although the two-state 
model cannot explain the molecular origins of such 
behavior, the fit of the data to eq (5) provides a 
basis of the intercomparison of temperature profiles 
in a methodical, reproducible, and quantitative 
manner. 

The values of the asymptotic slopes observed for 
the temperature profiles of intensity ratios cannot 
be interpreted on the basis of domains of small size 
alone. Even for the three conforming profiles in 
figure 5, the values of the slopes at high and low 
temperature are on the order of 0.5% to 1.2%, 
somewhat greater than the value of 0.46% ob- 
tained for the broadest domain size distribution re- 



1.60 r- 



1.40 



1.20 



1, 3/2-1 P 



D 2850/2885 
O 2940/2850 
a 2940/2885 



ported in table 1. The slopes are most readily inter- 
preted as the ability of the hydrocarbon chains of 
the lipid molecules to form gauche bonds near the 
acyl chain termini at the bilayer center without 
triggering the phase transition and, for the liquid 
crystalline phase, as a continuing disordering fol- 
lowing completion of the phase transition. 

Figure 6 compares the temperature profiles 
derived from the 7(2850 cm _1 )//(2885 cm" 1 ), 
7(2940 cm-y7(2850 cm" 1 ) and 7(2940 cm-')/ 
7(2885 cm -1 ) peak height intensity ratios for 1,3/2- 
lP-dipalmitoylcyclopentane-1,2,3- triol phosphat- 
ide acid, an analog of DPPC in which the glycerol 
backbone is replaced by a more rigid cyclopentane 
triol [13]. All three profiles give the same value for 
the phase transition temperature within the limits 
of the standard deviations obtained from the least 
squares fit. The values of AT, however, vary, and 
this variation can be used as a test of a structure- 
based model for the phase transition, since the 
7(2850 cm" 1 )/7(2885 cm" 1 ) ratio reflects primarily 
lateral chain-chain interactions, while the 7(2940 
cm _, )/7(2885 cm -1 ) ratio reflects trans-gauche 
conformational changes in addition to the lateral 
effects [14]. This suggests that different chain inter- 
actions as probed by different spectral regions may 
either exhibit different degrees of cooperativity in 
the phase transition, different energies for the tran- 
sition, or a combination of both effects. 

An even more striking comparison of tempera- 
ture profiles for the cyclopentanoid derivative of 
DPPC can be seen in figure 7. Simultaneous mea- 
surements of the C-H stretching and the C-C 
stretching mode regions showed a shift of 1.5° for 
the transition temperature for the C-C stretching 
profile, a shift which cannot be explained by the 
choice of intensities for forming the ratio as sug- 
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Figure 6-A comparison of the temper- 
ature profiles for three spectral fea- 
tures in the C-H stretching region 
for 1,3/2- lP-dipalmitoylcyclopen- 
tane-l,2,3-triol phosphatidic acid. 
For the 2850/2885 profile ( o). 
r =43.6 o C and Ar=0.6 "C; for 
the 2940/2850 profile ( o ), To 
= 43.4 °C and Ar=1.0°C; for 
the 2940/2885 profile ( a ), r 
=43.6 °C and Ar=0.8 °C. 
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Figure 7-A comparison of the temperature profiles for the C-C stretching 
and the C-H stretching mode regions for 1,3/2-lP-dipalmitoylcyclopen- 
tane-l,2,3-triol phosphatidic acid. The difference in T for the two spec- 
tral regions is 1.5 °C. Spectral data for both stretching mode regions were 
generated simultaneously. 
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gested by the discussion accompanying eq (6). Al- 
though the temperature intervals were too coarse 
to obtain precise values of A7\ the difference in T 
between the two profiles suggests that lattice ex- 
pansion precedes the intrachain trans-gauche iso- 
merization as the bilayer undergoes the phase 
transition. (The change in the intensity ratio in the 
C-H stretching region is not only a function of the 
degree of trans-gauche isomerization but also the 
interactions between neighboring hydrocarbon 
chains [1].) It is just this type of molecular-scale 
information that provides insight into the molecu- 
lar mechanism of the phase transition. Also, with- 
out a statistical estimate of the uncertainty in T , it 
would not be possible to assess whether the ob- 
served temperature difference is statistically signifi- 
cant. 

Finally, returning to figure 6, an asymmetry can 
be observed in the shapes of profiles which cannot 
be explained by a thermodynamic equilibrium 
model. The curvature at low temperature appears 
to be gradual, while at high temperature the curva- 
ture is quite sharp. Such an asymmetry is sugges- 
tive of a coalescing of domains of liquid crystalline 
lipid to give a narrower distribution of domains of 
larger size toward the end of the transition. This 
behavior is similar to that implied by the model of 
Freire and Biltonen [15], and observed by Marsh et 



al., using electron spin resonance techniques to ex- 
amine the DPPC phase transition [16]. 



Conclusions 

We have outlined a systematic approach for ex- 
tracting phase transition temperatures and van't 
Hoff enthalpies from temperature profiles derived 
from Raman spectral peak height intensity ratios. 
Within the degree of precision of the experimental 
data; namely, temperature intervals of 1-2 °C, tem- 
perature measurement accuracy of 0.1 °C, and peak 
height intensity measurement accuracy of 1 to 5%, 
the simple, equilibrium, two-state model implied by 
eq (5) fits most of the experimental data to within 
measurement error. Alternatively, with the above 
degree of precision of the experimental data, it will 
not be possible to distinguish between more de- 
tailed models. As a next step, therefore, in the sys- 
tematic evaluation of phase transition behavior 
through Raman spectroscopic studies, the preci- 
sion of the temperature and intensity measurements 
must be improved by an order of magnitude while 
decreasing the temperature intervals throughout 
the phase transition region by an order of magni- 
tude. 
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The advantage of deriving an analytical expres- 
sion for the Raman intensity behavior as a function 
of temperature lies in the development of statisti- 
cally based assessments of competing models. Data 
such as the exact profile curvature immediately 
preceding and following the phase transition and 
hysteresis effects can be used for model testing 
only when the data are reproducible and the mea- 
surement accuracy is controlled. A fruitful ap- 
proach is to combine improvements in the model 
with improvements in the experimental proce- 
dures, while using the statistical fit of analytical ex- 
pressions to the data as a basis for testing both the 
model and the adequacy of the data. 

We appreciate the limitations of the present 
model. Although the two-state model conveys no 
information concerning the molecular nature or 
origin of the phase transition, the vibrational spec- 
troscopic data provide a means for directly moni- 
toring structural changes within the lipid bilayer 
over a large temperature range. The spectroscopic 
data are thus complementary to the calorimetric 
data. Moreover, the extraction of thermodynamic 
information from the spectroscopic studies pro- 
vides confirmation that the spectroscopic and 
calorimetric studies are responding to the same 



phenomena. This aspect is particularly important in 
the study of lipid systems where sample history in- 
fluences the manifestation of conformational 
changes. As the next step in the development of the 
model, further detail concerning boundary effects 
between domains should be included along with a 
variation of the domain size distribution through 
the phase transition. 

With an increased number of parameters in an 
analytic expression for the phase transition profiles 
comes a demand for greater data density through- 
out the phase transition region. We have already 
indicated the need for decreased temperature inter- 
vals. The comparison of temperature profiles ob- 
tained from different spectral regions promises to 
provide extremely useful tests for molecular-scale 
descriptions of the phase transition and the atten- 
dant molecular dynamics. 



The course of the research described in this pa- 
per was strongly influenced by discussions with 
Dr. James S. Vincent and Dr. Terrell L. Hill. Their 
assistance is gratefully acknowledged. 



APPENDIX A 



Development of an Analytical Expression 
for Describing the Effect of Temperature on 
the Spectrum of a Molecular System in 
Terms of a Two-State Model. 



The total Raman scattered intensity 7 at a partic- 
ular frequency v can be expressed as: 



I v =l Av N A +I Bv N B t 



where I Av and Is* represent the scattering intensi- 
ties per molecule of molecules in each of two states 
A and B respectively and N A and N B are the num- 
ber of corresponding molecules per scattering vol- 
ume. The presence of other substances which may 
contribute to the overall scattering intensity will 
not affect the derivation that follows, since an ex- 
pression of the form 



L=IavNa +I&NB +I Cv N', 



(A2) 



where N' is the number of other scattering 
molecules, can be rewritten as 



(Al) I v =I Av N A +IB.NS +I Cv N'(N A +N B )/N 



=Ia*N a +U*Nb- 



(A3) 



I'a v and 7^,, as defined, are functions, therefore, of 
the composition of the sample as well as the intrin- 
sic scattering properties of the molecules of interest 
in states A and B. The temperature dependence of 
I v , however, will be determined by the temperature 
dependence of N A and N B . The purpose of the fol- 
lowing mathematical treatment is to determine an 
explicit analytical expression for this temperature 
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dependence. And since the treatment deals only 
with the numbers N A and N B , it is also applicable to 
absorption spectra in concentration ranges where 
the transmitted intensity or its logarithm is a linear 
function of concentration; that is, where Beer's 
Law is obeyed. 

We let A^ be the number of domains of size n of 
which N^ are in state A and N„ B are in state B. The 
number of molecules in states A and B are then 



N A =2nN n . 



and 



N B = 1nN nB . 



(A4a) 



(A4b) 



In a closed, one-component system in contact 
with a heat bath at pressure constant pressure p, 
the differential expression for the internal energy of 
the system is 

dE = TdS-pdV+2 iL„N n dn+2 jx„dN n (A5a) 

n n 

where 



li„=(\/N n )(dE/dn) s , KK 
and 

jx n =(dE/dN„) SfKn ^ m 



(A5b) 



(A5c) 



and where m refers to domains of all sizes other 
than n. The terms \i„N n dn and \i„dft n represent a 
separation of the total chemical potential of the sys- 
tem into two factors: one factor associated with 
changing the number of molecules in a given do- 
main and the other with changing the number of 
domains of a particular fixed size. In considering 
the effects of temperature on the number of do- 
mains in state A or state B, processes involving d« 
will be excluded. That is, only processes involving 
&ft nA and dN nB , where dA > ^ 4 = — &ft nB , will be al- 
lowed. 

The Gibbs free energy G for a system consisting 
of a mixture of domains of varying but distinguish- 
able size in one of two states, A or B, can then be 
written as: 

G = 2 £jV^ + /i fl A* (A6) 

n 

where 

fr=lN n =I.N nA +ft nB (A7) 

n n 

is the total number of domains in the system and 



where the chemical potentials fi^ and fi nB are the 
free energies of domains of size n in states A or B 
respectively. The condition for equilibrium be- 
tween states A and B is then: 



AG =2 £„ B -p^=0. 



(A8) 



We define the partition coefficients for the system 
of domains as: 



Z — H Z^Zrf 
n 

where 



__C 



M 



Z„,=-* 



(N n y. 



(A9) 



(A10) 



and £ ni is the partition function for a domain of 
size n in the state /. The chemical potential can also 
be defined in terms of the Gibbs free energy as: 



ix ni =(dG/d^ ni )T,p. n ,N m 



(All) 



and the Gibbs free energy can be expressed in 
terms of the partition function Z as: 



G = -kT\nZ+pV. 



(A12) 



Combining equations (A8) through (A 12) and in- 
voking Stirling's approximation for the logarithm 
of factorials gives: 



(A13) 



N nB 


-Aco/*r 


where 




AG n ° = - 


-kT\Jf- 

inA 



(A14) 



is the standard state free energy change for do- 
mains of size n. The number of molecules in state A 
is then given by: 



N A =Znti nA =Z 



nN n 



n 1+e 



-AG"/*r 



and in state B by: 

N B = 2nft nB =:2 



1+e + AGB/*Z- 



(A 15a) 



(A 15b) 



The temperature dependence of the intensity of a 
Raman spectral feature, which depends on N A and 
N B , is contained in the terms &G°/kT. As defined, 
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AG° is the difference in free energy between two 
domains of size n , one in state A and one in state B. 
For the simplest approximation to this difference, 
we let AG^ =n AG° A and AG Q nB =n AG° B where AGS 
and AGjj are the standard state free energies of indi- 
vidual molecules in states A and B respectively. 
This approximation is not exact because changes in 
the free energy of the domains will be affected by 
the domain boundaries. That is, we have defined 
the free energy of domains strictly for the internal 
degree of freedon associated with the order-disor- 
der phase transition of the lipid hydrocarbon 
chains. The free energy changes associated with 
this phase transition involve not only intramolecu- 
lar but intermolecular interactions as well. Individ- 
ual molecules at the boundaries between domains 
in different states will possess different intermolec- 
ular interaction energies in comparison with 
molecules within the domains themselves. The fac- 
torization of the partition function Z in eq (A9) into 
a product of single domain partition functions is 
also one that ignores boundary effects by assuming 
that the energy of interaction between domains is 
negligible. AG° and AG* as defined represent an 
average value of the free energy per lipid molecule 
for an average domain size. The definition of what 
constitutes an individual molecular species is flex- 
ible. A molecular species could refer to a complete 
lipid molecule, a single hydrocarbon chain or even 
a single -CH 2 -CH2-CH 2 -triad, the smallest unit for 
defining a trans or gauche conformation. The im- 
portant point is to ensure that the single unit free 
energy, AG° or AGs, is consistent with the defini- 
tion of that single molecular unit. In the remainder 
of this discussion, the unit will be taken as a single 
lipid molecule consisting of two chains. 

Substitution of AG n ° A =«AG£ and AG„° B =/iAGg 
into eq (A 15a) yields: 



N A =I 



nft n 



1+e 



-nAG»//?r- 



(A16) 



Note that a molar scale has been adopted for the 
exponential term so that AG is the free energy 
change per mole of lipid molecules, R is the gas 
constant, but n retains its definition as the number 
of molecules per domain. If n is large and if AG 
differs from 0, the exponential value is eq (A 16) 
rapidly approaches or infinity depending on the 
sign of AG and exhibits the discontinuity charac- 
teristic of a simple, first order phase transition. For 
smaller values of n (typically < 100), the system 
exhibits a more continuous transition from state A 
to state B and can be described in a manner 
analogous to a two-component system in chemical 



equilibrium. Since we know from experiment that 
pure lipid bilayers exhibit sharp phase transitions, 
we anticipate that n will be sufficiently large to 
allow expansion of AG°/T in a power series about 
T=T the temperature at which N A =N B : 



+ B(AG°/T)/dT 
J —l 



AG°/T=(AG°/T) 

+ 1/2 d\AG°/T)/dT 



\fl- T !°' 



JT-T a f. (A17) 

1 —Jo 



We define the standard state for each domain as the 
one for which N A =N B , so that AG°=0 at T=T Q . 
The derivative of AG°/T with respect to T can be 
evaluated from the Gibbs-Helmholz equation ex- 
pressed in terms of differences in partial molar en- 
thalpies; namely, 



d(AG°/T)/dT= -AH/T 1 



(A18) 



where AH=H B —H A . From this expression the sec- 
ond derivative of AG°/T with respect to T can be 
calculated: 

d 2 (AG°/T)/dT 2 =2AH/T 3 -(l/T 2 )(dAH/dT) 

(A19) 

Identifying dH/dT as the molar specific heat at 
constant pressure, C p , substituting eqs (A 18) and 
(A 19) into eq (A 17), and rearranging, the Taylor 
series expansion of AG°/T becomes: 

AG /T=-(&H/To)(T-T )/T 

+(AH/To-AC p /2)[(T-T )/To] 2 (A20) 

In general, phase transitions for lipids have been 
observed to be quite sharp with complete conver- 
sion occurring over a range of less than 10 K. 
Thus, at phase transition temperatures on the order 
of 300 K, (T-To)/T <.03 for r-r <10 K. Un- 
less AC P is an order of magnitude larger than AH/ 
T , the second term of eq (A20) can be ignored. 

Wilkinson and Nagle have measured the heat ca- 
pacities of several lipids within 5 K above and be- 
low their main phase transitions [17]. As the phase 
transition temperature is approached from lower 
temperatures, C p increases dramatically over the in- 
crease in C p observed for normal alkanes [18]. 
Above To, C p is higher for the lipid systems than 
for the alkanes and decreases with increasing tem- 
perature. However, AC, does not appear to in- 
crease significantly and may actually be lower than 
the value determined for the solid-liquid phase 
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transitions in alkanes. The anomalously high values 
of C p for lipids can be expected to occur when the 
phase transition is abnormally sharp; namely, over 
a temperature range of a degree or two where 
(T— T )/T <.005. Consequently, we assume that 
the temperature dependence of AG°/T can be ade- 
quately described by the first term in eq (A20). 

Substitution of the first term of eq (A20) into eq 
(A 16) gives for the number of molecules in state A , 



N A =2nft n /(\+e ,n ) 



(A21) 



where, for simplicity, we have defined a reduced 
temperature t as 



t = (T- T Q )/D = MI(T- T )/R Tl 
and the scale factor D as 
D=RT%/bH 

Similarly, 

N B = ZnN n /(l+e-' n ) 



{All) 



(A23) 



(A24) 



The magnitude of D may be estimated by letting 
A#^10 kcal/mole and 7^330 K so that D=*10 K. 
For n = 100 and T- T > 1 K, N A <N/\00. 

Equations (A21) and (A24) are still too complex 
to be used for a least squares analysis of the temper- 
ature dependence of spectral intensities in the 
vicinity of a phase transition. The next step in the 
development is to replace the summations in eqs 
(A21) and (A24) with some suitable average over 
the distribution of domain sizes. The distribution of 
domain sizes is given by N„ as a function of n, 
where N n satisfies the following conditions: 



n 

and 



(A25) 



(A26) 



and where the moments of the distribution are de- 
fined as 

<n k >=(\/ft)2ft n n k (A27) 

n 

The substitution of <n > into eq (A21), 

N A =fr<n >/(! +e' <n> ), (A28) 



would imply that <n > k = <n k > for all A: , which 
would be inaccurate for all but the narrowest of 
distributions. A better approximation is 



Ar<=JV<«>/(l+e/ 



</! 2 >/<n> 



(A29) 



Equations (A28) and (A29) both give the correct 
asymptotic behavior when compared with eq 
(A21). Both expressions have inflection points at 
t=0. However, eq (A29) also gives the correct 
slope at t=0, namely — A r <« 2 >/4, whereas eq 
(A28) gives -JV<iV> 2 . 

Since eq (A29) is congruent with eq (A21) at 
/ =0 and ± oo and since neither eq (A29) nor (A21) 
has any minima or maxima other than -f oo and 
-co, we conclude that eq (A21) may be approxi- 
mated by a single term, namely 



N A =N/(\+e' n '«) 



(A30) 



where N^N<n> and /t efr =<« 2 >/</i >. Simi- 
larly, eq (A24) can be approximated by 



N B =N/(l+e-"'<«) 



(A31) 



where N and n e{t have the same values as in eq 
(A30) in order to conserve the total number of 
molecules in the system. In practice, N and /i efr will 
be adjusted to give the best overall agreement be- 
tween eqs (A30,A31) and (A21,A24). 

The accuracy of this single-term approximation 
can be tested by assuming a distribution function 
for n and comparing eq (A30) directly with eq 
(A21). This was done for a set of lognormal distri- 
butions as described in the main body of this paper. 

Equations (A30) and (A31) can now be used in 
eq (Al) to describe the intensity of a Raman spec- 
tral feature as a function of temperature: 



1 = 



Ia v 



:+ 



(l+e') r (l+e-') 



(A32) 



where, as before, t=(T— T )/D, but where the ex- 
plicit reference to /i e(T has been absorbed into the 
definition of D, namely, 



D = {Rn)/LHn t „ 



(A33) 



and where explicit reference to the number density 
TV has been absorbed into the definitions of the scat- 
tering intensities I Av and 7^,. 
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APPENDIX B 



The Least Squares Fit. 



A computer program has been written to analyze 
the temperature dependence of Raman spectral 
data by fitting the data by the method of least 
squares to eq (5). Several features of this program 
are worth noting. Equation (5) is a linear function 
of the parameters A, B, a and b and a nonlinear 
function of T and D. Thus, the parameters of Eq 
(5) are evaluated using an iterative approach in 
which r is expanded in a Taylor series for each of 
the parameters about its currently estimated value. 
The series is terminated at its linear terms. The val- 
ues of rob,— ^0 are used as the dependent variables 
and the derivatives of r with respect to each of the 
parameters as the independent variables. From a 
linear fit of r obs — r^ to the derivatives of r, the 
corrections to the parameters are obtained and the 
process is repeated using the new values of the 
parameters until the changes in the parameters for 
each iteration are significantly lower than the stan- 
dard deviations of the parameters. 

The least squares fit is a weighted fit. If we as- 
sume that errors in the measured intensities are ran- 
dom, uncorrelated and are drawn from a single, 
normally distributed population of errors, and if we 
let s be the estimate of the standard deviation of 
that population, then the variance of a particular 
value of an intensity ratio is given by 



s\r n )=(dr n /dl x fs\l x )+(dr n /dhYs\h) 



=s 2 (\+rh)/ll 



where 



r\i—I\/h and s 2 (/,)=s 2 . Thus, the errors in the in- 
tensity ratios are not normally distributed but in- 
stead, depend on the magnitude of r and the 
intensity in the ratio denominator. Correspond- 
ingly, each ratio is given a weighting factor in the 
least squares fit proportional to the inverse square 
of its uncertainty, namely, 



*F=/f/(l+r 2 2 ). 



(B2) 



viation returned by the fit should be a measure of 
the uncertainty of the Raman intensities. Compari- 
son of the standard deviation of the fit with esti- 
mates of the uncertainty in the measured intensities 
demonstrated that eq (5) is capable of fitting the 
experimental data to within currently obtained 
measurement error. 

The algorithm for the least squares analysis fol- 
lows a Gram-Schmidt orthonormalization proce- 
dure. This procedure has been extensively tested 
[19] and a comparison between calculations using 
single and double precision gave identical results 
(to the degree of precision reported) indicating that 
the computation is free from numerical artifacts. 

The major problem encountered in the analysis 
of the temperature profiles is associated with those 
situations in which the observed data are insuffi- 
cient to determine the apparent phase transition 
temperature T and the scaling parameter D. The 
remainder of this section will deal with those situa- 
tions. 

The parameters T and D are determined solely 
from those data near T— T , whereas A,a,B and b 
are determined primarily by data from the high and 
low temperature asymptotes of eq (5). The range 
around T for which D and T can be determined is 
the range in which r deviates significantly from the 
asymptotes A +a(T—T ) and B+b(T— T ). A de- 
viation is defined as "significant" by the number of 
standard deviations by which a measured ratio dif- 
fers from each of the asymptotes. In particular, for 
a measured ratio to contribute to the determination 
of D and T , we require 

|r,— [a +a(T i -T )^>fs/Vw i (B3a) 



and 



(Bl) ^-(B+biTi-To^fs/Vw, 



(B3b) 



By setting the proportionality factor to unity rather 
than 1/s 2 as suggested by eq (Bl), the standard de- 



where r, is the value of the intensity ratio at the 
observed temperature T h A, a, B, b and T are the 
values of the parameters determined by the least 
squares fit; W t is the weight of the /th datum used 
in the least squares fit; s is the standard deviation 
returned by the least squares fit; and / is a multi- 
plicative factor to allow for a more or less stringent 
test. The quantity, s/VfVj, is the estimated uncer- 
tainty in the experimental value of the /' ,h ratio. 
Equations (B3a,b) state that for the purpose of con- 
tributing to the evaluation of T and D the value of 
an observed ratio must differ from both asymptotic 
limits by an amount greater than some multiple/of 
the expected measurement uncertainty at that tem- 
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perature. If the conditions of eqs (B3a,b) are met, 
we say that a measurement falls within the transi- 
tion interval. Equations (B3a,b) define the tempera- 
ture interval; the number of experimental data 
falling within this interval determines the reliability 
with which D and T can be estimated. In particu- 
lar, the number of data falling within this interval is 
the appropriate number to use in calculating the 
number of degrees of freedom for the estimation of 
confidence intervals from a student's / distribution. 

If no measured temperature falls within the tran- 
sition interval, neither T nor D can be estimated by 
a least squares fit of the data to eq (5). If one mea- 
sured value falls within the interval, only T can be 
estimated, but not D. Furthermore, the value of the 
uncertainty in T returned by the least squares fit 
must be multiplied by a factor of 12.7, the value of 
t from the student's / distribution for one degree of 
freedom and confidence limits of 95%. If two mea- 
sured values fall within the interval, both T and D 
can be estimated, but the estimates of their uncer- 
tainties must be multiplied by 4.3, the value of / 
from the student's t distribution appropriate for 
two degrees of freedom. Should more than two 
measured values fall within the transition interval, 
T Q and D can be determined with increasing reli- 
ability. 

When less than two points fall within the transi- 
tion interval, bounds can be placed on the values of 
T and D, even though T and D cannot be deter- 
mined precisely. In order to place bounds on T 
and D, bounds must be placed on the transition 
interval. The transition interval can be estimated 
by combining eqs (B3a,b) with eq (5) to give, ig- 
noring terms in a and b, 

T X =F X D, 



where 



F x =2ln( 



B-A 



Vw 



fi 



-i) 



(B4) 



and Wis, the average weight of a measured value in 
the vicinity of the phase transition temperature. 
The transition interval T x is the temperature inter- 
val over which r varies from A +a(T— T )±fs/ 
Vw to B+b(T-T )+fs/Vw, where the upper 
sign is chosen if A <B and the lower sign if B <A. 
If no measured temperature falls within the tran- 
sition interval, as tested by eqs (B3a,b), we know 
that 

0<T x <T d (B5) 

where T d is the temperature interval between the 



successive data pair above and below the phase 
transition temperature. Thus, from eq (B4) 



0<D<T d /F x 



(B6) 



In calculating the value of A , a, B and b, D is held 
fixed at its average value in the range defined by eq 
(B6); namely 



D = T d /2F x , 



(B7) 



and T is given as the midpoint between the two 
measured temperatures spanning the phase transi- 
tion interval. Depending on the precision of the 
Raman intensity measurements relative to the dif- 
ference between A and B, F x typically varies from 
2 to 9. Most of the time a value of/=2 is used in 
eqs (B3a,b) and (B4). Substitution of &T=4D into 
eq (B6) yields 



bT=*RTl/(tJln tn )<AT d /F x 



and 



n e „>F x Rn/(LHT d ). 



(B8a) 



(B8b) 



If one measured temperature T t falls within the 
transition interval, 



0<T x <2T d , 



(B9) 



where 2T d = Tj+\ — T t — 1 is twice the average 
data interval of the three measured temperatures 
spanning the transition interval. As with the case of 
no measured temperature falling within the inter- 
val, D is assigned the average value implied by the 
range of eq (B9); namely, 



D = T d /F x , 
and AT is reported as 
±T<%T d /F x . 



(BIO) 



(BH) 



T is estimated by the least squares fit of all the data 
to eq (5) but is determined, almost exclusively, by 
the value of r at T t relative to A and B. 
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